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TEN  ORBIT  DETERMINATIONS  FOR  COSMOS  472  ROCKET 

by 

Doreen  M.  C.  Walker 


SOttlARY 

e  472  rocket,  1972-04B,  has  been  determined  from  461 
observations  at  tan  epochs  during  the  last  two  months  of  its  life.  The  first 
eight  orbits  were  satisfactory,  with  average  standard  deviations  equivalent  to 
220  m  in  cross-track  distance  and  75  m  in  radial  distance. 

The  variation  of  perigee  height  has  been  analysed  to  obtain  four  values  of 
density  scale  height.  The  first  three  showed  good  agreement  with  CIRA  1972. 

Unfortunately  the  geometry  of  the  orbit  was  such  that  the  decrease  in 
inclination  did  not  yield  an  exact  valua  of  atmospheric  rotation,  contrary  to 
expectations. 
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UNLIMITED 
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INTRODUCTION 


Cosmos  472  rocket,  I972-04B,  vu  Che  rocket  which  fired  Che  satellite  Coeaos  472 
into  orbit  on  1972  January  25.47.  After  separation,  1972-04B  had  an  orbital  period1  of 
102  minutes,  perigee  height  193  km,  apogee  height  1526  km  and  inclination  82.0°.  The 
rocket  remained  in  orbit  for  101.58  days,  decaying  on  1972  May  6.05. 

With  this  initial  orbit,  1972-04B  promised  to  be  useful  for  determining  the 

atmospheric  rotation  rate  from  the  decrease  in  orbital  inclination.  The  orbit  of 

2 

19 72-0 4A  had  been  used  previously  to  obtain  a  value  of  atmospheric  rotation  rate,  A  , 
with  values  of  inclination  from  OS  Navy  orbits;  but  the  value  of  A  obtained  was  not 
considered  to  be  of  a  very  high  standard. 


The  orbit  of  1972-04B  has  been  determined  from  461  observations  with  the  aid  of  the 
BAE  orbit  refinement  program^  PROP  6,  at  ten  epochs  during  the  last  two  months  of  its  life 
Unfortunately  however,  the  first  six  values  of  inclination  obtained  were  at  a  time  when 
(■>  ,  the  argument  of  perigee,  was  near  270°  and  consequently  there  was  little  change  in 
inclination,  so  chat  A  could  not  be  accurately  determined.  As  w  moved  towards  180° 
and  a  change  in  inclination  became  apparent,  the  orbit  also  passed  through  15th-order 
resonance.  As  the  satellite  was  rapidly  decaying  by  this  time,  it  was  not  possible  to 
calculate  the  perturbation  due  to  resonance  from  the  sparse  data  available  and  also 
accurately  evaluate  a  value  of  A  . 

The  main  reason  for  determining  the  orbit  having  proved  unfruitful,  the  perigee 
height  was  snalysed  to  determine  the  atmospheric  density  scale  height.  Four  values  of 
scale  height  were  determined  and  the  results  are  presented  in  this  Memorandum. 

2  THE  0BSESVAXI0MS  AMP  ORBITS 

The  orbit  of  1972-04B  has  been  determined  at  ten  epochs  using  517  observations. 

The  type  and  nuaber  of  observations  used  on  each  orbit  are  given  in  Table  1.  The 

Table  I 


Sources  of  observations  used  on  each  orbit 


Orbit 

No. 

Source  of  observations 

Total 

Visual 

us 

Navy 

Finland 

British 

radar 

1 

53 

20 

8 

81 

2 

25 

50 

12 

8 

95 

3 

5 

32 

37 

4 

7 

54 

61 

5 

41 

41 

6 

46 

46 

7 

10 

32 

8 

so 

8 

1 

19 

25 

45 

9 

4 

24 

28 

10 

5 

20 

9 

34 

Total 

no 

338 

21 

49 

518 

observations  fall  into  four  groups;  the  visual  Observations  made  by  volunteer  observers 
now  reporting  to  the  Barth  Satellite  Issearch  Unit  at  die  University  of  Aston,  the 


US  Havy  observations  supplied  by  tha  US  Haval  Research  Laboratory,  cha  Finnish  observa- 
tions  nada  by  tha  thaodolita  at  Jokioinan,  Finland,  and  observations  nada  by  British 
radars. 

Tha  orbital  alaaants  at  aach  of  tha  tan  apochs  ara  listad  in  Tab la  2  on  paga  5  with 
tha  standard  dsviations  balov  aach  valua;  tha  apoeh  for  aach  orbit  is  at  00  hours  on  tha 
day  indicatad.  Tha  PROP  6  nodal  fits  tha  naan  anomaly  M  by  a  polynomial  of  tha  form 

M  -  Mq  *  M,t  ♦  M2t2  ♦  Mjt3  ♦  M4t4  ♦  Mjt5  (I) 

vhara  t  is  tha  tima  naasurad  from  apoch  and  tha  nunfcer  of  M-coaf f iciants  usad  dapands 
on  tha  drag.  Bast  rasults  vara  obtainad  by  using  Mq  to  M^  ,  tha  full  complement  of 
coafficiants  in  tha  modal,  on  two  of  tha  tan  orbits,  two  raquirad  only  Mq  to  Mj  ,  and 
tha  ramaining  six  raquirad  Mq  to  M^  . 

Tha  orbits  fit  tha  obsarvations  in  a  satisfactory  nannar  on  tha  first  sight  orbits, 
with  e  ,  tha  paranatsr  indicating  tha  msasura  of  fit  ranging  batvaan  0.45  and  1.0.  Tha 
last  two  orbits  wars  nora  difficult,  as  tha  satallita  was  decaying  rapidly  and  there  ware 
lass  obsarvations  available. 

Tha  PROP  program  rejects  obsarvations  which  do  not  fit  wall,  and  461  of  the 
original  518  obsarvations  ware  accepted  in  tha  final  orbits.  Tha  residuals  of  tha 

A 

obsarvations  have  bean  obtainad  using  tha  ORBS  computer  program  and  sent  to  tha 
observers. 

The  accuracies  of  salactad  observing  stations,  with  five  or  more  obsarvations 
accepted  in  tha  orbit  determination,  ara  listad  in  Table  3  along  with  tha  number  of 
accepted  and  rejected  obsarvations.  Three  of  tha  visual  stations  have  had  their  rms 
residuals  revised  because  during  tha  orbit  determination  a  few  of  their  observations  had 
thair  standard  deviations  increased  in  order  to  keep  them  in,  rather  than  allowing  them 
to  be  rejected.  Tha  revised  sat,  omitting  these  obsarvations  (never  more  than  15X), 
gives  a  batter  impression  of  tha  observers'  accuracy. 

3  INCLIMAIIOH 

For  a  high-drag  satallita,  atmospheric  rotation  is  tha  most  important  force  per¬ 
turbing  tha  inclination  i  ,  so  if  tha  change  in  i  is  msasured  accurately,  and  other 
perturbations  ara  removed,  cha  rotation  rata  of  tha  upper  atmosphere  in  tha  region  near 
tha  satellite's  perigee  can  be  determined. 

Tha  values  of  inclination  from  Table  2  are  plotted  in  Fig  1  after  removal  of  tha 
zonal  harmonic,  Jj  j  ***  luniaolar  perturbations.  Tha  tonal  harmonic  and  luni solar 
perturbations  ware  removed  by  using  the  PROD*  computer  program  with  1-day  integration 
steps,  and  the  Jj  2  tes serai  harmonic  perturbation  by  using  tha  valua  recorded  on 
the  PROF  printout.  Perturbations  due  to  earth  and  ocean  tides  were  expected  to  be 
smeller  than  tha  standard  deviations  of  tha  values  of  i  ,  so  were  ignored. 

The  theoretical  variation  of  i  due  to  atmospheric  rotation  on  be  calculated  for 
a  series  of  values  of  atmospheric  rotation  rate,  A  ,  using  the  computer  program  MZfflf 


7 


based  on  equation  (32)  of  Ref  6.  Unfortunately  the  value  of  argument  of  perigee  m 
palm  through  270°  just  after  orbit  3,  vith  the  result  thst  the  dsersas*  in  i  over  the 
first  six  values  is  very  snail  and  there  is  no  chance  of  determining  an  accurate  value  of 
A  .  Over  the  last  month,  however,  the  inclination  falls  by  0.025°;  but  by  a  further  mis¬ 
fortune  the  satellite  passes  through  15th- order  and  31:2  resonances  during  this  time. 

The  theoretical  curve  vith  A  •  0.7  gives  the  best  fit  to  the  values  of  inclination  in 
Fig  1. 

This  value  of  A  is  unexpectedly  low  and  so  the  effect  of  a  meridional  wind  was 
investigated.  The  ROT ATM  program  also  allows  for  the  effect  of  meridional  winds  but  for 
a  satellite  with  an  inclination  of  82°  the  effect  is  small,  being  <0.001°  for 
u  ■  0.1  rev/day  where  v  is  the  south-to-north  rotation  rate.  It  would  seem  therefore 
that  the  values  of  i  being  fitted  may  have  been  perturbed  by  one  or. both  of  the 
resonances,  and,  as  the  data  at  the  times  of  the  resonances  is  so  sparse,  the  perturba¬ 
tions  cannot  be  reliably  calculated. 

The  height  for  which  the  value  of  A  applies  is  |H  above  perigee,  where  H  is  the 
density  scale  height.  This  gives  a  value  of  220  km,  and  from  a  recent  review  of  upper 
atmosphere  winds ^ ,  the  lowest  value  of  A  is  0.78  at  this  height;  and  that  is  only 
obtained  when  the  satellite  is  experiencing  morning  and  susmer  conditions.  This 
satellite,  1972-04B,  was  not  in  sumner  conditions,  the  last  four  orbits  being  determined 
between  15  April  and  4  May,  so  at  best  could  only  be  described  as  having  a  tendency  to 
suamer.  The  local  time  over  this  final  stage  varied  from  14  h  to  10  h,  so  two  of  the 
orbits  were  determined  in  morning  conditions.  From  Ref  7  a  low  value  of  A  (near  0.9) 
would  be  expected  for  a  satellite  with  the  local  time  and  seasonal  conditions  of  I972-04B, 
but  not  as  low  as  0.7. 

Another  explanation  could  be  that  before  the  last  epoch  the  satellite  was  perturbed 
in  its  orbit  by  the  escape  of  some  residual  propellant,  thus  causing  an  increase  in  the 
value  of  inclination.  There  is  some  support  for  this  theory  in  the  analysis  of  the 
perigee  height  in  section  5. 

4  PERIGEE  HEIGHT 

The  values  of  perigee  distance  a(l  -  e)  are  plotted  as  circles  in  Fig  2,  using 
the  values  of  a  and  a  from  Table  2.  If  the  odd  zonal  harmonic  and  lunisolar  pertur¬ 
bations  are  removed  from  a(l  -  e)  ,  using  PR0D^,  the  remaining  variation  should  be  a 
steady  decrease  as  a  result  of  air  drag  alone.  These  corrected  values  of  perigee 
distance,  Q  ,  for  the  ten  epochs  of  Table  2  are  plotted  as  crosses  in  Fig  2,  with  a 
smooth  curve  drawn  through  the  crosses. 

The  actual  perigee  height,  yp  ,  above  the  Earth's  surface  is  found  by  first 
restoring  the  zonal  harmonic  and  lunisolar  perturbations  to  the  smoothed  values  of  Q  , 
giving  smoothed  values  of  perigee  distance,  a(l  -  Smoothed  '  ctMQ  Cram  these  values 
subtracting  the  local  Barth's  radius,  Rp  ,  at  latitude  4p  ,  given  by 

R_  •  R  -  21.38  sin2  * 

P  P 
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I 
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* 


i 


•t 


vh«r%  R  it  tht  Earth's  tquttoritl  radius,  6378.14  lea;  and  finally  adding  tha  small 

g 

amount  by  which  the  actual  path  of  tha  satallita  departs  from  an  exact  allipsa  ,  dr  , 

P 

whara 

dr  -  1 

P 


.388 


sin 


■ 


Thus,  sines 


sin2  4 


.  2  .  .2 

sin  i  sin  u 


7p  “  4(1  "  *> smoothed  ‘  6376*75  +  (20‘96  *  l¥|)  «in2  -  • 

Tha  values  of  y^  calculatad  at  tha  tan  epochs,  using  tha  above  aquation,  are  plotted  in 
Fig  3  and  a  s Booth  curve  has  bean  drawn  through  tha  points. 

5  TENSITY  SCALE  HEIGHT 

Tha  perigee  distance  cleared  of  zonal  harmonic  and  1 unisolar  perturbations,  as 
given  by  Q  ,  gradually  decreases  under  the  influence  of  air  drag,  and  the  decrease  is 
proportional  to  tha  density  scale  height,  H  .  The  values  of  Q  ,  given  by  the  crosses 
in  Fig  2,  show  this  steady  decrease  and  values  of  H  can  be  found  from  the  changes 
in  Q  . 

o 

The  theoretical  equation  for  the  variation  of  Q  is 


2*' 


4ae  ~  a 


sin  i  cos  2w  ♦  0 


(l  H  B2  e'2\l 


(3) 


where  H  is  the  density  scale  height,  is  the  value  of  B  at  perigee,  Bj  is  the 
value  of  B  at  a  height  l.SB^  above  perigee,  and  e'  is  the  allipticity  of  the  atmos¬ 
phere  ,  taken  as  the  Barth's  allipticity,  0.00333.  This  equation  for  Q  is  valid  for 

aa/H  >  3  . 

P 

Average  values  of  Q  were  calculated  over  a  time  interval,  At  ,  long  enough  to 
ensure  that  an  accurately  measurable  change,  AQ  ,  in  Q  has  occurred.  The  values  of 
AQ/At  serve  as  values  of  Q  in  equation  (3),  and  averaged  valuas  of  Mj,  M2,  e,  a,  i 
and  cos  2w  over  the  corresponding  time  At  were  used  to  calculate  values  of  Bj  from 
equation  (3) . 

Four  values  of  B}  were  calculated  and  they  are  plotted  as  circles  in  Fig  4,  and 
the  time  over  which  they  are  averaged  is  indicated.  The  corresponding  average  values  of 
height  y( (-  yp  ♦  1.58^)  are  given  at  the  top  of  Fig  4. 

The  values  of  B|  plotted  as  crosses  in  Fig  4  are  the  valuaa  of  density  scale 
height  obtained  from  the  COSPAR  International  Bafarence  Atmosphere  1972  (CIRA  1972)  for 
heights  y|  and  appropriate  exospheric  temperature. 

The  agreement  between  the  first  three  values  of  B(  evaluated  here  and  the  CIRA 
values  is  very  good,  as  the  CIRA  values  would  be  expected  to  have  an  error  of  up  to  I0X. 
The  last  value  of  B{  determined  from  the  orbital  data  appears  to  be  too  low.  Mo 
geophysical  reason  can  be  found  for  this,  as  agreement  with  the  CIRA  value  would  require 


1 1 1 


perigee  diatance  to  b«  5  kn  lover.  Tbit  difference  u«u  to  bo  for  too  lagge  to  bo 
to  orbital  errora ,  which  art  of  order  300  a,  ao  again  tha  poaaibility  that  aoata 
irnt  propellant  eacaped  fxoa  tha  rocket,  thua  changing  ita  orbit,  auat  bo  given  aa  tb 
noat  likoly  explanation. 

6  COMGLB8IOBS 

Tan  orb ita  of  1972-048  have  boon  computed,  eight  being  of  aatisfactory  accuracy. 
They  fail  to  give  a  good  value  of  A  ,  aa  had  boon  hopod,  becauao  the  valuaa  of  e  were 
unfavourable;  but  throe  good  aatinatOa  of  atmoapheric  acalo  height  have  boon  derived. 


REFERENCES 


Author 

D.G.  King-Hele 
J.A.  Pilkington 
H.  Hiller 
D.M.C.  Walker 

D.G.  King-Hele 
D.M.C.  Welker 

R.H.  Gooding 

D.W.  Scott 

G.E.  Cook 


D.M.C.  Walker 

D.G.  King-Hele 
D.M.C.  Walker 

T.  Koxai 

D.G.  King-Hele 
D.M.C.  Walker 


AH,' 

a'  ■  i;  • 

OH  TQ  uOM 


Title,  etc 

The  RAE  table  of  Earth  satellites  1957-1980. 

Macmillan  Press,  London  (1981) 

Upper-atmosphere  zonal  winds:  variation  with  height  and  local 
time. 

Planet.  Space  Sri.,  25,  313-336  (1977) 

RAE  Technical  Report  76055  (1976) 

The  evolution  of  the  PROP  6  orbit  determination  program,  and 
related  topics. 

RAE  Technical  Report  74164  (1974) 

ORES:  a  computer  program  for  the  analysis  of  residuals  from  PROP. 
RAE  Technical  Report  69163  (1969) 

PROD,  a  computer  program  for  predicting  the  development  of  drag- 
free  satellites.  Part  1:  theory. 

RAE  Technical  Report  71007  (1971) 

[Celestial  Meohanios,  301-314  (1973)] 

Cosmos  462  (1971-106A):  orbit  determination  and  analysis. 

Phil.  Trans.  Roy.  Soo.  A,  292,  473-512  (1979) 

RAE  Technical  Report  78089  (1978) 

Upper-ats»sphere  zonal  winds  from  satellite  orbit  analysis. 

RAE  Technical  Report  (to  be  issued) 

The  motion  of  a  close  Earth  satellite. 

Astronam.  Joum.,  64,  367-377  (1959) 

Air  density  at  heights  near  150  km  in  1970,  from  the  orbit  of 
Cosmos  316  (1969-108A). 

Planet.  Space  Sri.,  _19,  1637-1651  (1971) 

RAE  Technical  Report  71129  (1971) 

CIRA  1972  (Cospar  International  Reference  Atmosphere  1972). 
Akademie-Verlag,  Berlin  (1972) 


"TO  ARC  NOT  rr^r ? ' 

'  ?  v - ■  •  v  ■  i  ■  f 

LHCIAL  OHOANIRAiioNP 


Local  time  at  perigee 


04  24  20 


81.97 


81.96  |- 

Inctinat  ion 
deg 

81.95 


81.94 


81.93 


15th-order 

resonance 


31:2 

resonance 


1972 

Mar  14  Mar  24  Apr  3  Apr  13  Aor  23  M 


Date  -  modified  Julian  day 


Values  of  inclination  cleared  of  perturbations 
with  theoretical  curve  for  A*0.7 


Date 'modified  Julian  day 


Fig  2  Values  of  a(1-e) 


(1440 


